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Abstract. The problem of the limit shape of large alternating sign matrices 
(ASMs) is addressed by studying the emptiness formation probability (EFP) 
in the domain-wall six-vertex model. Assuming that the limit shape arises in 
correspondence to the 'condensation' of almost all solutions of the saddle-point 
equations for certain multiple integral representation for EFP, the limit shape 
of large ASMs is found. The case of 3-enumerated ASMs is also considered. 



1. Introduction 

An alternating sign matrix (ASM) is a matrix of l's, O's and — l's such that in 
each row and column, all nonzero entries alternate in sign, and the first and the 
last nonzero entries are f 's. In a weighted enumeration, or q-enumeration, ASMs 
are counted with a weight q k , where k is the total number of —l's in each matrix. 
There are many nice results concerning ASMs, mainly devoted to their various 
enumerations; for a review see, e.g., book [1]. 

In this paper, we address the problem of the limit shape of large ASMs. The 
problem comes out from the fact that in ASMs their corner regions mostly contain 
O's while in the interior there are many nonzero entries. As the size of ASMs in- 
creases, the probability of finding entries equal to 1 (or —1) in the corner regions 
vanishes, while in the central region such probability remains finite. When consid- 
ering very large ASMs in an appropriate scaling limit, (e.g., when large matrices 
are scaled to a unit square) , it is common knowledge that such regions are sharply 
separated by a smooth curve. 

In essence, such curve is an Arctic curve for ASMs, similar to the Arctic Circle 
of domino tilings of large Aztec diamonds [2] . Since the Arctic curve determines the 
shape of the internal, or "temperate" , region of ASMs, we call it here the limit shape 
of ASMs (cf. [3]). More generally, one can address this problem for q-enumerated 
ASMs. The case q = 2 corresponds to the Arctic Circle of the domino tilings. Here 
we obtain explicit expressions for the limit shapes of 1- and 3-enumerated ASMs. 

To treat the problem, we exploit the one-to-one correspondence between ASMs 
and configurations of the six-vertex model, which has been found and efficiently 
used in papers [4-6] . This correspondence takes place when the six- vertex model is 
considered with the so-called domain wall boundary conditions (DWBC). The six- 
vertex model with these boundary conditions has been originally introduced and 
studied in [7-9] . 



2 



F. COLOMO AND A.G. PRONKO 



a 




a 




c 

-1 



Figure 1. The six allowed arrow configurations, their weights, 
and the corresponding entries of ASMs. 
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Figure 2. A possible configuration of the six- vertex model with 
DWBC and the corresponding ASM, for N = 5. 



Using the standard description of local states in terms of arrows (see, e.g., [10]), 
DWBC mean that the model is considered on a square lattice of N vertical and N 
horizontal lines, where the arrows on all external horizontal edges point outward, 
while the arrows on all external vertical edges point inward. Then, the 0's in ASM's 
entries correspond to the vertices of weights a and 6, while both l's or — l's entries 
correspond to vertices of weight c, as shown on Figure 1. Figure 2 shows a possible 
configuration of the six-vertex model DWBC and the corresponding ASM. 

To count ASMs, the weights a, 6, c are to be put all to the same value, e.g., 
a = b = c = 1. In ^-enumeration, ASMs are counted by taking a = b = 1 and 
c = y / q 1 since in each configuration the vertices of type five and six come in pairs, 
in addition to N vertices of type six being always present due to DWBC. 

The main tool which allows us to study limit shapes of the domain wall six- 
vertex model is a particular, non-local, correlation function, the so-called emptiness 
formation probability (EFP). In paper [11] some multiple integral representation 
for EFP has been derived. Here, we the study the multiple integral representation 
for large size of the lattice (corresponding to large dimension of ASMs). Assuming 
that the limit shape occurs in correspondence to the peculiar situation where almost 
all roots of the saddle-point equations condense to the same value, we derive the 
Arctic curves in the two cases of values of weights matching 1- and 3-enumeration 
of ASMs. This 'condensation hypothesis' holds rigorously in the free-fermion case, 
considered in [12], where it allowed us to give an alternative derivation of the Arctic 
Circle (the limit shape of 2-enumerated ASMs). 

Our paper is organized as follows. In the next Section we recall results of [11] 
on multiple integral representations. In Section 3 the condensation hypothesis is 
discussed, and a procedure for the derivation of limit shapes is explained. In Section 
4, we derive the limit shape of 1- and 3-enumerated ASMs. 
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2. Emptiness formation probability 

As in [11,12], we call emptiness formation probability (EFP) and denote by 
F^' s \ where r,s = 1,...,N, the probability of having all arrows on the first s 
horizontal edges from the top of the lattice, located between r-th and (r + l)-th 
vertical lines (counted from the right), to be all pointing left. 

With this definition, EFP measures the probability that all vertices in the top- 
left (N — r) x s sublattice have the same configuration of arrows, namely, with 
all arrows pointing left or downwards, or, equivalently, that these vertices are of 
type two (see Figure 1). This follows from the peculiarity of both the domain- 
wall boundary conditions and the six-vertex model rule of two incoming and two 
outgoing arrows at each lattice vertex (known also as the 'ice-rule'). 

In the language of ASMs the indicated property implies that all entries in the 
top-left (N — r) x s block are O's. The limit shape then arises in some appropriate 
scaling limit, and corresponds to some curve where EFP jumps from one to zero as 
the size of this (N-r)xs block increases. It is to be mentioned also that the study 
of the behaviour of EFP provides only the top-left 'quarter' of the ASM's limit 
shape; the remaining three quarters follow from simple symmetry considerations. 

In [11] several equivalent representations for EFP were given. For what follows 
we shall need two representations in terms of multiple integrals. To recall these 
formulae, we introduce some objects first. 

An important role in our considerations is played by the function hpf(z) = 
hpf(z; A, t), where A and t are to be regarded as parameters, 

. a 2 + b 2 -c 2 b 

A:= — , *:=-, 

2ab a 

while z is to be treated as a variable. This function, defined as a generating function, 
is a polynomial of degree (N — 1) in z, 

N 

h N (z):=J2 H N^ r ~\ Ml) = l- 

r=l 

Here the quantity = (A, t) is a boundary correlation function, introduced 
in [13]. Namely, it is the probability that the sole vertex of type six (having weight 
c), residing in the first row, appears at r-th position from the right. At t = 1 
the function Iin(z) has a special meaning as the generating function for refined 
^-enumerations of ASMs, with q and A related by A = 1 — q/2. Our derivation of 
the limit shapes below involves essentially the known explicit expressions for this 
generating function at some particular values of q. 
For s = 1, . . . , N, we define functions 

h NtS {zu...,z s )= H (zj-zk)- 1 det [z^izj - l) 5 - fe /ijv_fe + i(^)] . 

-*■ ■*■ 1<?'.K<S J 

l<j<k<s - 

These functions can be regarded as multi- variable generalizations of Iin(z), in the 
sense that Hn,i{z) = Iin{z). It can be easily checked that 

ftjv,s+i(zi) • • • > Zs, 1) = h,N, s (zi, ■ ■ ■ i z a ). (1) 

One also has 

h NjS+ i(zi, ...,z s ,0) = h N (0)h N -i,s(zi, ■■■,z s ). (2) 
Properties (1) and (2) are used in what follows. 
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We are now ready to turn to the multiple integral representations. In [11], the 
following multiple integral representation has been obtained 



x n t 2 ^2tAz J+ i ^ (zi --' z - )dzi --- dz - (3) 

Here Co denotes some simple anticlockwise oriented contour surrounding the point 
z = and no other singularity of the integrand. Formula (3) has been derived in 
[11], using the quantum inverse scattering method [14,15] and the results of papers 
[8,9,13]. This formula has also been discussed in [12] (see equation (3.6) therein). 
In [11], the following equivalent representation has been also given: 



(_l)»(«+i)/2 Za r r ' [{f _ 2tA)zj + 1] 



j?(r,s) _ K-i-) _ fs / i TT 

N ~ .c\C?.^\s„s(s-l)rS T„ "T„ 11 



ls-1 



S !(2 7 ri)«a«(«- 1 )c« J Co J Co ^ z v( Zj - 1)- 

x n f2 fc _ Laz, + 1 n ^-^ 2 

j.k=l ■> ■> l<j<k<s 

x h N . s ( Zl , . . . , z s )h SyS (u 1 , . . . , u s ) dz! ■ ■ ■ dz s . (4) 

Here Z s denotes the partition function of the six-vertex model with domain wall 
boundary conditions on s x s lattice, and 

Zi — 1 



*3 ■- 



(t 2 - 2tA)zj + 1 ' 



(5) 



For later use we emphasise that the integrand of (4) is just the symmetrized version 
of the integrand of (3), with respect to permutations of the integration variables 
z\,...,z s . This follows through the symmetrization procedure explained in [16], 
and some additional identity proven in [11]. 

We are interested in the behaviour of EFP in the so-called scaling limit, that is 
in the limit where r, s and N are all large, with the ratios r/N and s/N kept finite 
(and smaller than 1). By applying standard arguments of saddle-point analysis to 
representation (4), we obtain the following system of coupled saddle-point equations 



s(t 2 - 2tA) 



1 (t 2 — 2tA) Z j + 1 Zj Zj — Zk 



fc=i 



^ / 2At- t 2 z k t 2 z k \ 8 



fc=i 



t 2 ZjZk — 2AtZj + 1 t 2 ZjZk — 2Atzk + 1/ dzj 
t 2 - 2At + 1 d 



[{t 2 - 2At) Zj + l] 2 d Uj 



\nh StS (ui, ...,u s ) = 0. (6) 



In deriving these equations we have used the fact that quantities like ln/ijv iS are of 
order s 2 , and that their derivatives with respect to Zj's are of order s; all sub-leading 
contributions (estimated as o(s)) are neglected. 
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3. Condensation hypothesis and 'reduced' saddle-point equation 

To address the problem of limit shapes, one need, in principle, to be able to 
describe solutions of saddle point equations (6). This task however is a formidable 
one, at least because the last two terms in (6) are rather implicit in general. Never- 
theless, it turns out possible to outflank the obstacle provided that some hypothesis 
holds true. 

In [12], it has been shown that in the case of A = the limit shape corresponds 
to the situation where almost all roots of (6), as s — > oo, condense at point z = 1. 
More precisely, the condensation has been assumed first, and next, after obtain- 
ing the limit shape, the validity of this assumption has been confirmed. Here we 
conjecture that the same correspondence of limit shape and condensation of roots 
can be put forward also in a more general situation, since the main ingredients 
appear not to be specific to the A = case. The only difference is that now such 
assumption cannot be verified a posteriori (at least by simple methods), so we call 
it here 'condensation hypothesis'. 

When assuming condensation, in [12] we have relied on two observations which 
can be easily extended beyond the free fermion case. The first one, based on the 
physical interpretation of EFP as a correlation function which tests ferroelectric or- 
der macroscopically, is that in the scaling limit EFP has a step-function behaviour. 
The second observation is that some multiple integral, with just the same inte- 
grand as that appearing in representation for EFP, but with integration contours 
surrounding the point z = 1, is equal to one exactly. 

Namely, basing on (3), let us consider, for r, s = 1, . . . , N, the integral 



(r , s) _ (-iy r r ^ [(t 2 - it a) Zj + 1]^ 



x n t 2 *i aL+i h N ,M,...,z s )d Z1 ...dz s . 



l<j<k<s 

Here C-f is a closed contour in the complex plane enclosing point z = 1 and no 
other singularity of the integrand; the minus in the notation indicates negative 
(clockwise) orientation. We have 

IP S) = 1. (7) 



Indeed, performing integration in variable z s , and taking into account (1), identity 
(7) follows immediately by induction. Note that relation (7) implies the same 
relation for the integral with the integrand of (4), which is just the symmetrized 
version of the integrand of (3). 

Thus, we can explore again, just as in [12], the condensation of saddle-point 
solutions at the value z = 1, as a possible mechanism for the emergence of limit 
shapes. For this purpose we introduce here the concept of the 'reduced' saddle 
point equation, which allows one to treat the consequences of such condensation 
rather directly, and in the general setting of equations (6). Below we rely on 
results of papers [17,18], where the scenario of condensation of roots of saddle- 
point equations has been described in the context of random matrix models with 
logarithmic potentials (Penncr models, [19]). 
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To consider the scaling limit we define x and y such that 

N ~ r s , n „ 

x:=— — , y-=N' x,ye[0,l\. 

By 'condensation' we mean that, as s — > oo, almost all roots Zj of (6) have the value 
z = 1, in the sense that the amount of non-condensed roots is of order o{s). This 
implies that at s = oo the density of roots is just the Dirac delta-function S(z — 1), 
or that the Green function is G(z) = (z — 1) . To implement this condensation in 
(6), we pick up j-th equation and just set z^ = I for k ^ j. Writing simply z for 
this Zj, we arrive at the following equation 

y 1 - x t 2 y . v 1 9 

r n ; 7+ hm — — mh,N(z) = 0. (8) 

z - 1 z i 2 z - 2tA + 1 w^oo dz w w 

In deriving this equation from (6), we have used properties (1) and (2), in particular, 
from (5) we have Uk = as Zk = 1, and hence property (2) shows that the last term 
in (6) becomes just the derivative of a constant, thus vanishing. 

We call equation (8) the reduced saddle-point equation. The solutions of this 
equation give the non-condensed roots of (6) at condensation. As explained in [18], 
a necessary condition for condensation of roots of saddle-point equations like (6), 
is the presence of two coinciding non-condensed roots. In our case these roots must 
lie on the real axis, in the interval [l,oo). 

Summarizing, we arrive at the following recipe for the derivation of limit shapes: 
we require (8) to have two coinciding roots, which moreover must run over the 
interval [1, oo). Denoting the value of these two roots by w, we shall see below that 
as lo runs over the real axis from point z = 1 to point z — > oo, it parameterizes the 
top-left portion of the limit shape, from the top contact point (corresponding to 
u> = 1) to the left contact point (corresponding to u = oo). 

An ingredient which is necessary for this programme, is the knowledge of the 
last term in (8). Fortunately, function Iin(z) is explicitly known in some interesting 
cases. These are, for t = 1, the case of A = 1/2, corresponding to usual enumeration 
of ASMs, and the case of A = —1/2, corresponding to 3-enumerated ASMs; for 
generic t the function Hm{z) is also known for A = 0, corresponding to the six-vertex 
model in the frcc-fcrmion case (specializing further t = 1, one gets 2-cnumcratcd 
ASMs). 

To illustrate how the recipe is working, let us consider the case of A = 0. We 
set also t = 1 for simplicity, so that function Hm{z) — hjy(z; A, t) in this case is just 
/i/v(z; 0, 1) = [(z + 1)/2] N ~ 1 (see, e.g., [12]). The reduced saddle-point equation 
reads 

V x 1 — y 
— ^— -- + = 0. 

z — 1 z z + 1 

Denoting by g(z) the function in the left-hand side, we require g(z) = (z — ui) 2 g(z) 
where g(u) ^ 0. This can be implemented by solving the system of two equations 
g(u) = and g'(co) = 0, where the prime denotes derivative, for unknowns x and 
y. This gives 

Eliminating uj, we also have the equation for the limit shape 

4x(l - x) + 4y(l - y) = 1. (9) 
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Here x and y take values in interval [0,1/2]; equation (9) describes the top-left 
portion of the Arctic Circle, which is the limit shape of 2-enumerated ASMs. 

4. Limit shapes of 1- and 3-enumerated ASMs 

We start with considering the case of 1-enumerated ASMs. In this case we have 
A = 1/2 and t — 1, and the function Kn{z\ 5, l) is given by the formula (see, e.g., 
[20,21]) 



M«l.i)-.«(-™H 



We write here the hypergeometric function in such a way that the third parameter 
is positive, and larger than the second one, so that the Euler integral representation 
for Gauss hypergeometric function can be used to study the large N limit. 
Explicitly, the Euler integral representation gives the following expression 

M*! h 1) - |^ jf [-(1 - r)(l r + ZT)}"- 1 dr. 

The large N behaviour of this integral can be found via the standard saddle-point 
analysis, and we find that, as N — > 00, 

\nh N (z;\,l) = N\n[Av(l-v){l-v + zv)] + o{N), 

where 

2 - z - Vz 2 - z + 1 

v := . 

3(1 -z) 

The reduced saddle-point equation in this case reads 

' l-x + y + l-y/#-z + \ = ^ (10) 



z — 1 z z(l — z) 

The requirement that equation (10) has two coinciding roots gives the following 
parametric solution for the limit shape of large ASMs: 

2uj - 1 oj + 1 r 

x = 1 = , y = 1 = , u> € 1 , 00 ) . 

Eliminating u from this parametric solution we obtain that the limit shape is de- 
scribed by the equation 

4r(l - x) + 4y(l - y) + Axy = 1, x,ye[0,\]. (11) 

Formula (11) is the central result of our paper. 

The most refined numerical simulations of large ASMs (up to N = 1300) have 
been performed recently by Wieland [22], who also provided pictures comparing 
these numerical data with equation (11). As the size of ASMs increases, convergence 
to this curve is observed. 

Let us now consider the case of 3-enumerated ASMs (A = —1/2 and t = 1). 
The function Iin(z; — |; l) has been obtained in [20] (see also [21]). The following 
formulae are valid 

h , _ 1 . 1 v = Ul/2)(z+l)B 2m (z) for N = 2m + 2, 

N[Z ' 2 ' ' \(l/9)(2z + l)(z + 2)B 2m (z) for N = 2m + 3, 
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Figure 3. The limit shapes of g-enumerated ASMs for q = 1, 2, 3, 
given by equations (11), (9), and (14), respectively. 



where Bi m {z) is the following polynomial of degree 2m: 



B 2m {z) = 



m + 1 



3 m - 1 (2m + 3) 
m 



z m (z + 2) m 2 F 1 



— m, m + 2 
2m + 4 



z 2 -l 



5m— 1 



(2m + 3) 



z m (z + 2) 1 "- 1 2 Fi 



z(z + 2) 
m + 1, m + 2 
2m + 4 



z(z + 2) 



Again, in comparison with [21], we have written here the hypcrgeometric polyno- 
mials in such a way that the Euler integral formula can be directly applied. 

Evaluating the integrals through the saddle-point method, we obtain that, as 
N -> oo, 



\nh N (z;-\,l) = Nln 

The reduced saddle-point equation reads: 
V l-x y 



2(2z + l)(z + 2) 



9(^+1) 



+ o(N). 



z-l 



+ 



2z 2 + 4z + 3 



2 + z (l + z)(2 + z)(l + 2z) 



0. 



(12) 



The requirement that equation (12) has two coinciding roots gives us the fol- 
lowing parametric solution for the limit shape of 3-enumerated ASMs: 



7uj 4 + 14a; 3 + 19^ 2 + 12^ + 2 
(lu 2 + 2)(2u + 1) 2 (lu + 1) 2 ' 
(w - l) 2 (6w 4 + 16w 3 + 19w 2 + 16w + 6) 
3(w 2 + 2)(2lV + l)2( w + l)2 ' 



w e [l, oo). 



(13) 



Equivalcntly, one can search for the equation connecting x and y. In this way we 
obtain that the limit shape of 3-enumerated ASMs is described by the following 
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sextic equation: 

324 x 6 + 1620 x 5 y + 3429 x 4 y 2 + 4254 x 3 y 3 + 3429 x 2 y 4 + 1620 xy 5 + 324 y 6 

- 972 x 5 - 1458 x 4 y - 2970 x 3 y 2 - 2970 x 2 y 3 - 1458 xy 4 - 972 y 5 

- 6147 x 4 - 9150 a; 3 y - 17462 x 2 y 2 - 9150 xy 3 - 6147 j/ 4 

+ 13914 x 3 + 24086 a; 2 y + 24086 xy 2 + 13914 y 3 

- 11511a; 2 - 17258a;?/- 11511 y 2 

+ 4392 x + 4392 y - 648 = 0, a;, y e [0, ±] . (14) 

One can verify directly that this equation is indeed satisfied by x and y given by 
(13). 

Figure 3 shows plots of the limit shapes for 1-, 2- and 3-enumerated ASMs. 
The area of the temperate region decreases as q increases, in agreement with both 
analytical and numerical considerations [23-25]. 
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